Electron-phonon coupling in semimetals in a high magnetic field 
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We consider the effect of electron-phonon coupling in semimetals in high magnetic fields, with 
regard to elastic modes that can lead to a redistribution of carriers between pockets. We show that 
in a clean three dimensional system, at each Landau level crossing, this leads to a discontinuity in the 
magnetostriction, and a divergent contribution to the elastic modulus. We estimate the magnitude 
of this effect in the group V semimetal Bismuth. 



In a quantizing magnetic field, the electronic density 
of states is split into Landau levels, and in a clean three- 
dimensional system, there is a square-root singularity in 
the density of states as a function of energy. This is the 
origin of a host of magneto-oscillation effects, that in- 
ter alia provide a detailed tomography of fermi surfaces 
in metals. There are effects which are peculiar to multi- 
band rather than single-band systems. In a system with a 
single species of carrier, the chemical potential oscillates 
with field in order to keep the carrier number fixed. But 
with multiple pockets of carriers, there is an opportu- 
nity to redistribute carriers between pockets to minimize 
the energy with coupling generated cither by electron- 
electron or electron-phonon coupling. In this paper, we 
concentrate on the latter, and especially coupling to uni- 
form strain in the crystal. 

The study of large magnetostriction in semimetals goes 
back to Kapitzai and Shoenberg^, with work from the 
1960's onwards concentrating on quantizing fields^^— up 
to 20 T. Bismuth (Bi) has recently come to the fore in 
a series of experiments in higher field^i^— - the extreme 
quantum limit where only one or a few Landau levels are 
occupied - which have produced puzzling results that we 
will discuss in more detail below. Part of the purpose of 
this paper is to point out the extent to which electron- 
phonon coupling should be included in any analysis, and 
to detail which of the phenomena observed experimen- 
tally cannot be produced by magnetostrictive effects, in 
order that the arena to search for new physics is clear. 

We begin by a simple formulation of the problem and 
recapitulate general results. We then estimate the mag- 
nitude of the effects in Bi, and discuss this in the context 
of both older and more recent results. Since the defor- 
mation potential coupling at the band edge is not well 
established experimentally, we use ab-initio calculations 
of the band edges under strain to guide our estimates. 

For an illustrative discussion of the effects, wc will con- 
sider the case of two pockets because no particular extra 
physics is introduced by further generalization. The Lan- 
dau levels are then indexed by quantum number n and 
momentum k in the direction of magnetic field; again for 
simplicity we will assume the field to be oriented along 
a principal axis of the effective mass tensor. The carrier 



energies in the two pockets z = 1 , 2 are 
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where is the effective mass, ujic = eB/rriicC the cy- 
clotron frequency (see Fig. [T]). Holes are treated as hav- 
ing positive energies. A chemical potential is intro- 
duced to fix the carrier densities rii, and the band edges 

^ + li.aSa shift with Strain Sq,. While it is true 
that other parameters - for example the effective mass - 
arc strain-dependent, it is only the band-edge shift that 
contributes in a singular fashion, as seen below. 

We specify two cases of interest: 

1. Electron and hole pockets: ni = n2 constrains 
equal numbers in the two pockets. Here, — 
-f 7s adjusts the (equal) numbers of electrons 
and holes so that the strain is emptying electrons 
from one pocket to the other. 




FIG. 1. Sketch of the Landau levels of an electron and hole, 
where the chemical potential is indicated by a dashed line. 
The lower figure represents hole energies as positive, the con- 
vention used here. 
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2. Two degenerate electron pockets: ni + 712 = n, 
with the total carrier number fixed. Here, £1.2 ^ 
ei.2 ± 7s, so that the strain breaks the symmetry 
between the pockets. 

Of course, the type of strains in the two cases will be 
different. 

We now compute, conventionally, the total electronic 
energy and the total occupancy. In each Landau level, 
the density of states (per unit energy, per unit area) is 



where now 
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and, therefore, the carrier density and the total electronic 
energy in each pocket become 
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We define 



1 /'^ 



+ (n + 1/2) hhJic + ^iS , 



(5) 



i = {he/ eBY^'^ is the magnetic length, and Q{x) is the 
Heaviside function. 

It is convenient to rescale lengths in units of the mag- 
netic length £ and energies in units of the cyclotron en- 
ergy of one of the pockets tuoic- Thence, in dimensionless 
units, we have 
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The total energy of the system includes the strain energy 
and the electronic terms 



(8) 
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We are interested in following the total energy as a 
function of strain and magnetic field (which acts as a 
tuning parameter). To determine this, Eq. ([6]) fixes the 
chemical potential, fed in to determine the total energy 
via Eq. ([7]). In general this is not analytically tractable 
except in special cases (e.g. particle hole symmetry, when 
/^ = 0). _ 

However, it is clear that the contribution of strain to 
the electronic energy is generally small, and therefore it 
is appropriate to separate out the singular contributions 
that occur near where Landau levels cross through the 
chemical potential from a smooth background, viz. 
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So 



1 AS 
7~b" 



(10) 



is the tuning parameter (magnetic field), which van- 
ishes when the n*^ Landau level of interest empties: 
AS = Bn — B. (We are thus confined to Case I, for the 
moment.) The smoothly varying part of the energy can 
be incorporated into regular magnetostriction that will 
lead to a small spontaneous strain, that weakly renormal- 
izes the effective stiffness, Kcs, and shifts the parameter 
sq. We re-parameterize to take this smooth variation into 
account. (And parenthetically we note that other strain- 
dependent parameters, such as the effective mass, may 
be taken into account in the same way.) The last term 
is non-analytic, and needs to be treated separately from 
the background variation. 

By straightforward minimization, one finds 



for ?/ < , 
for ?7 > , 
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FIG. 2. (Color online) The effective elastic modulus from 
Eq. (|13[) . induced strain from Eq. (jlip . and Landau level oc- 
cupancy from Eq. (|15p . as a function of ry = AB/B. The 
latter is also proportional to the torque on the sample, when 
there is angular variation of Bn{9). 
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where we introduced the parameters 77 = A_B / B and 
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Corresponding to the kink in the spontaneous strain, 
there is a singular contribution to the effective stiffness, 
i.e. 
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After minimization, the internal energy becomes 
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and the Landau level occupancy is given by 



An = 2 — = ^ 
orj 7 
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For large strains, the above approaches the Landau level 
occupancy result obtained when the induced strain is ig- 
nored, i.e. 



An,=o = 2 
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where the K and 7 dependences cancel out. The size 
of the effects induced by strain is controlled by the small 
dimensionless parameter /3, which is (re-dimensionalizing 
the other parameters) 
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The range in reduced field where the elastic corrections 
are observed is of order Ary ~ 7/?/ (huJc)- 

The above results are depicted in Fig. [5] and represent 
the generic case of non-degenerate pockets (i.e. Case I). 
Case II, where there are degenerate pockets, is of interest 
in case the mechanism predicts an instability, breaking 
the symmetry between the pockets. By inspection, it is 
clear that it does not, because the energy increases as a 
power of the tuning parameter with an exponent greater 
than unity (3/2). Were the exponent less than unity, 
there would be a discontinuous jump in density. 

but in dimensions less than 2 an instability is expected. 

When the magnetic field is not aligned with a princi- 
pal axis of the pockets, there is a torque, r = dU/d9. 
At the Landau level edge this acquires a discontinuous 
contribution, i.e. 



dU d\ogB„ 



(18) 



Note that there is a contribution to the torque in the 
absence of coupling to strain, owing to the diamagnetic 
contribution of the carriers in the pockets. In the case 
considered here, the contribution at small excess fields 
arises because of the physical transfer of carriers from one 
pocket to another: the strain coupling reduces the torque 
and changes the behavior from square root to linear. 

We now turn to a discussion of elemental Bi. Bi 
is a classic semimetal, having a hole pocket around 
the L-point on the (111) surface of the Brillouin zone, 
and three electron pockets around the T-points on the 
(Til), (111), (111) directions. One may interpret the 
band structure by realizing that the rhombohedral struc- 
ture of Bi is derived from a distorted simple cubic struc- 
ture: there is a dimerization along the (111) direction, 
which forms a puckered sheet (though the dimerisation 
is small, with a difference of bond to back-bond lengths 
of 12% ) and opens small gaps at the fermi energy; this 
is accompanied by a rhombohedral distortion from the 
ideal angle of 60° to 57.2° which breaks the symmetry 
between the L and T points on the Brillouin zone sur- 
face, leading to band overlap and a semimetal. Conse- 
quently, the band edges of the conduction and valence 
points are expected to be very sensitive to rhombohe- 
dral strain oriented along (111); orthorhombic strain will 
break the symmetry of the electron pockets. 

Since the early work of Kapitza and Shoenberg, Bi 
has long been known to have very large magnetostric- 
tionii^, and the effects of quantizing magnetic fields have 
been studied by several authors^"—, with Michenaud et 
al^ presenting data up to 19 T and an analysis similar 
to the one above. 

Recently, a number of studies have explored high- 
field (> 20T) characteristics of Bi, though not explic- 
itly the magnetoelastic effects. Nernst and Hall mea- 
surementsi revealed unexpected structure attributed to 
the hole pocket at fields large enough that the chemi- 
cal potential is expected to be inside the lowest Landau 
level. Torque measurements as a function of angle (which 
are sensitive to the electron pockets) revealed sharp (and 
hysteretic) structure with strong angular dependence^; 
subsequently the position (field, angle) of these anomalies 
has been identified— with electron pocket Landau level 
structure. The physics of this anomaly remains unidenti- 
fied, but is consistent with breaking of the electron valley 
degeneracy. An alternative explanation— was proposed 
to explain the coincident observation of Hall plateaus at 
certain magic angles. 

We focus on estimating the order of magnitude of the 
coupling the soft "C44" rhombohedral shear, which is ap- 
proximately the longitudinal strain along the trigonal 
axis. We calculate the electron-phonon coupling from 
first principles by varying the rhombohedral angle and 
calculatingi^ the energies of the band edges at the L 
and T points, yielding 7 « 10 ± 1 eV. At 20 T, we 
estimate /3 « 10~^, consistent with the numerical esti- 
mates of Michenaud et al^ using parameters extracted 
by Walthcr— , and with their measurements at fields up to 
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19 T. Though the induced strains are small, they are sig- 
nificant. The range in reduced field where the corrections 
arc seen is P^/hu)c ~ 0.1 — 0.5, so that the conventional 
formulae for n{^), E{^) are invalid; the modulation in 
the pocket density is of order Kfi/^, about 10^^ cm^"^ 
at 20 T using the numbers above ~ around 1% of the 
zero-field carrier density. 

We thus remark that the magnetoelastic effects in Bi 
are expected to be unusually large, especially at high 
fields. As well as the predictions of "single-particle" 
physics, the sensitivity of strain to Landau- level occu- 
pancy would make such experiments useful to study 
interaction-driven transitions in this system, should they 
exist. Preliminary high- field measurements show indeed 
substantial magneto-elastic quantum oscillations at fields 
exceeding 20 T (see Ref. [l3l ). The importance of such 
measurements is that they are sensitive to the total elec- 
tron and hole numbers in the pockets, and are intrinsi- 
cally easier to interpret than transport data. We predict 
a singularity in the sound velocity at the Landau-level 
edge in a clean system. 

With the exception of the high field Ncrnst measure- 
mentsi on the fractionally filled hole pocket, there is as 
yet no strong evidence for correlation effects in Bi at 
other than close to integer Landau level fillings. Nonethe- 
less, there are several puzzles in the data that have been 
taken as evidence for the modulation of carrier densities 
in the valleys. It is important to distinguish between ef- 
fects that indicate a spontaneous instability of the system 
from those that can be derived mechanically. 

Consequently, the electron Landau-level lines that rise 
steeply at high fields, that were identified by the field 
H2{0) (Fig. 3 of Ref. in torque measurements, are not 
special for their existence (which follows the most rc- 
ceniiiSiii^ though not earlier calculations), but for their 
shape. The observed hysteresis near zero angle is not 
compatible with single-particle effects as discussed here, 
though the approximate shape of the anomaly is consis- 
tent with Fig. [2] for angles in the range of 2-4 degrees 
where no hysteresis is reported. 

The observed (sec Fig. 8 of Ref. [1) reversal of the 
sign of the torque at magic angles for rotations in the 
trigonal-bisectrix plane is consistent with the crossing of 
(different) Landau levels of the doubly-degenerate and 
singly-degenerate electron pockets, with the sign-change 
arising from the opposite angular-dependence of the Lan- 
dau level shifts. Associated with this effect, it was noted 
that the Hall resistance pxy - in this geometry believed 
to be dominated by the holes - is field-independent over 
a range of several Tesla. In contrast, for rotations in the 
trigonal-binary plane (when only a single electron pocket 
is involved), the Hall resistance shows a pronounced an- 



gular minimum but continues to grow with field. We note 
that the magnetoelastic modulation of the carrier density 
due to rhombohedral strain will modulate electron and 
hole densities equally; this is a mechanism for the elec- 
tron Landau-level crossings to moderate transport by the 
holes, aside from any direct contribution of the electron 
carriers to transport coefficients. Since An cx t, one may 
speculate that the near cancelation of the torque anoma- 
lies along the trigonal-bisectrix axis might be reflected in 
Pxy, in contrast to the single Landau level physics in the 
binary direction. 

Finally, we must make a distinction between the elas- 
tic strain induced effects discussed here (at long wave- 
lengths, (7 = 0) and the possibility of finite— g phonon 
softening and instabilities due to field-induced nesting, 
much studied in the context of field-induced charge and 
spin density waves^^. At zero magnetic field, the Bi elec- 
tron and hole pockets do not nest, and so a density- 
wave instability is suppressed. At large fields, the re- 
duced dimensionality of the dispersion will give rise to 
one-dimensional-like nesting features and hence an in- 
stability. Note, however, that (at least in weak coupling) 
the transition temperature will generically scale with a 
coupling constant A as Tc oc e~^/^, where A oc rriiz via 
its dependence on the density of states. The very light 
mass of carriers in Bi drives such an instability to very 
low temperature. 

In conclusion, we have made a quantitative estimate of 
magnetoelastic coupling effects in a degenerate scmimctal 
in the quantum limit, predicting the detailed shape of 
anomalies that should be measurable in Bi. In contrast to 
two dimensions, we find that electron-phonon effects per 
se do not produce instabilities, but lead to discontinuities 
in derivatives of the elastic strain and a singularity in the 
modulus. 

We believe that carrying out such measurements in Bi 
at high fields ~ and preliminary data has been obtainedi^ 
- would be very helpful in unraveling the complex physics 
that has emerged from other recent experiments^^— on 
this system. 
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